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Abstract. At variance with most semi-analytic models, in the Anti-hierarchical Baryon 
Collapse scenario (Granato et al. 2001, 2004) the main driver of the galaxy formation and 
evolution is not the merging sequence but are baryon processes. This approach emphasizes, 
still in the framework of the hierarchical clustering paradigm for dark matter halos, feedback 
processes from supernova explosions and from active nuclei, that tie together star formation 
in spheroidal galaxies and the growth of black holes at their centers. We review some recent 
results showing the remarkably successful predictive power of this scenario, which allows 
us to account for the evolution with cosmic time of a broad variety of properties of galaxies 
and active nuclei, which proved to be very challenging for competing models. 

Key words. Cosmology: theory - galaxies: formation - galaxies: evolution - quasars: gen- 
eral 



1. Introduction late, by the merger of disk/bulge systems made 

primarily of stars. 

The establishment in the 1990's of the hier- 
archical clustering paradigm led to the devel- On the other hand, it has long been known 
opment of various semi-analytic models for that stellar populations in elliptical galaxies 
galaxy formation sharing the basic assumption are old and essentially coeval (Sandage & 
that the main driver shaping the structure and Visvanathan 1978; Bernardi et al 1998; Trager 
morphology of galaxies is gravity. In this see- et al. 2000; Terlevich & Forbes 2002). A color- 
nario, the gas cools and form stars following magnitude relation is also well established: 
the collapse of dark matter halos. In a cold dark brighter spheroids are redder (Bower et al. 
matter cosmology, small objects form first and 1992). The widely accepted interpretation is 
merge together to make larger ones. This see- that brighter objects are richer in metals and 
nario then implies that large ellipticals form the spread of their star formation epochs is 
small enough to avoid smearing of their col- 
Send offprint requests to: G. De Zotti ors. The slope of this relation does not change 
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with redshift (Ellis et al. 1997; Kodama et al. 
1998) supporting this interpretation. The star 
formation history of spheroidal galaxies is mir- 
rored in the Fundamental Plane (Djorgovski & 
Davies 1987; Dressier et al. 1987) and in its 
evolution with redshift. Elliptical galaxies ad- 
here to this plane with a surprisingly low or- 
thogonal scatter (~ 15%), as expected for a ho- 
mogeneous family of galaxies. Recent studies 
(e.g. Treu et al. 2002; van der Wei et al. 2004; 
Holden et al. 2004, 2005) suggest that ellipti- 
cals, both in the field and in clusters, follow 
this fundamental relation up to z i 1, consis- 
tent with the hypothesis that massive spheroids 
are old and quiescent. 

Although these data have mostly to do with 
ages of stars and may leave open the issue of 
the epoch at which ellipticals were assembled 
to their present size and mass, they motivated 
our group to take a different view, i.e. to in- 
vestigate the possibility that gas processes have 
a key role in driving the formation and evo- 
lution of large galaxies, still in the framework 
of hierarchical clustering of dark matter halos 
(Monaco et al. 2000; Granato et al. 2001). This 
new approach emphasized the role of feed- 
backs, first of all from supernovae, which re- 
lease large amounts of mechanical energy, ca- 
pable of unbinding the gas in weakly bound, 
low mass systems, but also from the active nu- 
clei which were found to be ubiquitous in the 
centers of spheroidal galaxies (Kormendy & 
Richstone 1995; Magorrian et al. 1998). These 
feedbacks can actually reverse the order of for- 
mation of visible galaxies compared to that of 
dark halos: large galaxies form their stars first, 
while the star formation is delayed in smaller 
halos (Anti-hierarchical Baryon Collapse, or 
ABC, scenario; Granato et al. 2001, 2004). At 
the same time, processes associated to the star 
formation activity have a profound effect on 
the evolutionary history of nuclear activity. 

In this paper we will briefly review some 
recent results obtained in the framework of the 
ABC scenario, focussing in particular on pre- 
dictions for the evolution of nuclear activity. 
We adopt a spatially flat cold dark matter cos- 
mology with cosmological constant, consis- 
tent with the Wilkinson Microwave Anisotro- 
py Probe (WMAP) data (Bennett et al. 2003): 



Q. m = 0.3, n b = 0.047, and Q A = 0.70, 
H Q = 70km s" 1 Mpc~', <x 8 = 0.84, and an in- 
dex n — 1 .0 for the power spectrum of primor- 
dial density fluctuations. 

2. Correlations between SMBHs and 
their host galaxies 

In the last several years various empirical re- 
lationships between super-massive black hole 
(SMBH) masses and different properties of 
their host galaxies have been derived (but 
see Novak et al. 2006). These include cor- 
relations with the stellar velocity dispersion 
(Ferrarese & Merritt 2000; Gebhardt et al. 
2000; Tremaine et al. 2002; Ferrarese & Ford 
2005), with the mass in stars (Haring & Rix 
2004), with the bulge luminosity (McLure & 
Dunlop 2002; Marconi & Hunt 2003; Bettoni 
et al. 2003), with light concentration (Graham 
et al. 2001), and with the dark halo mass 
(Ferrarese 2002; Baes et al. 2003). In all cases, 
BH masses are related to the (generally old) 
bulge stellar population, not to the younger 
disk (Kormendy & Gebhardt 2001; Kormendy 
& Ho 2000; Salucci et al. 2000). 

Shankar et al. (2006) exploited updated es- 
timates of the total and baryonic mass func- 
tions of galaxies, of their luminosity function, 
of their velocity dispersion function, and of the 
mass function of central black holes to derive 
relationships between halo masses and stellar 
masses, stellar velocity dispersions, and black 
hole masses (Mbh)- Particularly relevant, in the 
present context are the relationship between 
Mh and M star , both in solar units: 

M star (M^xlO 11 ) 31 
1O 1O M ~ ' l+(M h /3xlO u ) 2 - 2 ' 1 ' 

and between Mbh and M/,: 

M BH (MJ2.2X10 11 ) 3 - 95 
10 7 M o 1 +(M A /2.2x lO 11 ) 27 ' 

3. Velocity Dispersion Function and 
Virial Velocity Function 

As first pointed out by Loeb & Peebles 
(2003), a comparison of the local Velocity 
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Dispersion Function (VDF) with the Virial 
Velocity Function (VVF) can provide interest- 
ing hints on the structure formation process. 
Accurate determinations of the VDF of early- 
type galaxies have been obtained by Sheth et 
al. (2003) and Shankar et al. (2004), based on 
a large sample (~ 9000 E/S0 galaxies) drawn 
from the SDSS (Bernardi et al. 2003). As for 
dark matter halos, it is convenient to define a 
'virial' velocity, equal to the circular velocity 
at the virial radius (Navarro et al. 1997, NFW; 
Bullock etal. 2001): 



Since, given the virialization redshift, W de- 
pends only on M vir (V vir oc A/V 3 ) the VVF 
can be straightforwardly derived from the mass 
distribution function of spheroidal galaxies, in- 
tegrated over the virialization redshifts. 

Cirasuolo et al. (2005), assuming that all 
massive halos (2.5 x 1O H M < M vir < 2 x 
10 13 M o ) virializing at z > 1.5 yield spheroidal 
galaxies or bulges of later type galaxies, 
showed that the VVF accurately matches the 
VDF for a constant ratio of the velocity dis- 
persion cr to the virial velocity, cr/y vir 0.55, 
close to the value expected at virialization if it 
typically occurred at z i 3. This is a remark- 
able result, since the VVF depends only on the 
evolution of dark matter halos, while the VDF 
is affected by the physics of baryons. 

A substantial stability of the halo circu- 
lar velocity after the fast accretion phase was 
found by Zhao et al. (2003) in high-resolution 
N-body simulations, even though the halo mass 
increases by a substantial factor. But the central 
velocity dispersion may be affected by merg- 
ing events and dissipative baryon loading. This 
strongly suggests that dissipative processes and 
later merging events had little impact on the 
matter density profile, consistent with the dy- 
namical attractor hypothesis (Loeb & Peebles 
2003; Gao et al. 2004). It may also suggest 
that, for massive galaxies, the mass assembly 
is largely complete at substantial redshifts, as 
also indicated by the Spitzer Space Telescope 
observations by Papovich et al. (2006). Direct 
evidence that the galaxy stellar mass function 



does not evolve significantly since z - 12 has 
been presented by Bundy et al. (2006). 

These results may be relevant also for 
the interpretation of the M^-cr correlation, 
which probably originated at substantial red- 
shifts, when the old stellar population formed. 
If so, both M B h and cr should not have changed 
much since then. 



4. Black hole growth 

Another important ingredient of the evolution- 
ary history of quasars is how the active nuclei 
acquired their mass. Although radiative accre- 
tion must play a role (Soltan 1982), it is not 
necessarily the main mechanism for the BH 
mass growth. In the hierarchical assembly sce- 
nario for massive galaxies a buildup of SMBHs 
by coalescence of BHs associated to merging 
subunits is naturally expected. 

On the other hand, the analysis by Shankar 
et al. (2004; see also Marconi et al. 2004), 
exploiting up-to-date luminosity functions of 
hard X-ray and optically selected AGNs, has 
shown that the local SMBH mass function is 
fully accounted for by mass accreted by X- 
ray selected AGNs, with bolometric correc- 
tions indicated by current observations and 
a standard mass-to-light conversion efficiency 
e ^ 10%. An unlikely fine tuning of the pa- 
rameters would be required to account for the 
local SMBH mass function accomodating a 
dominant contribution from 'dark' BH growth 
(due, e.g., to BH coalescence). It may be noted 
that the local SMBH mass function can be 
rather accurately assessed. Estimates from ei- 
ther the velocity dispersion function of galax- 
ies coupled with the M B h-o" relation or from 
the galaxy luminosity function coupled with 
the MsH-isph relation of McLure & Dunlop 
(2002) agree very well (cf. Fig. 1 1 of Shankar 
et al. 2004). 

The visibility time, during which AGNs are 
luminous enough to be detected by the cur- 
rently available X-ray surveys, is found to be 
of 0.1-0.3 Gyr for present day BH masses 

M BH - 106 " 1O 9 M . 
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5. The Granato et al. (2004) model 

5.1. Overview of the model 

While referring to the Granato et al. (2004) pa- 
per for a full account of the model assump- 
tions and their physical justification, we pro- 
vide here, for the reader's convenience, a brief 
summary of its main features. 

The model follows with simple, physically 
grounded recipes and a semi-analytic tech- 
nique the evolution of the baryonic compo- 
nent of proto-spheroidal galaxies within mas- 
sive dark matter (DM) halos forming at the 
rate predicted by the standard hierarchical clus- 
tering scenario within a ACDM cosmology. 
The main novelty with respect to other semi- 
analytic models is the central role attributed 
to the mutual feedback between star forma- 
tion and growth of a super massive black hole 
(SMBH) in the galaxy center. 

The idea that SN and QSO feedback play 
an important role in the evolution of spheroidal 
galaxies has been pointed out by several au- 
thors (Dekel & Silk 1986; White & Frenk 
1991; Haehnelt, Natarajan & Rees 1998; Silk 
& Rees 1998; Fabian 1999). Granato et al. 
(2004) worked out, for the first time, the sym- 
biotic evolution of the host galaxy and the cen- 
tral black-hole (BH), including the feedback. 
In this model, the formation rate of massive 
halos (2.5 x 10 n M o < M vir <2x 10 13 M o ) is 
approximated by the positive part of the time 
derivative of the halo mass function (Press & 
Schechter 1974, revised by Sheth & Tormen 
2002). The gas, heated at virial temperature 
and moderately clumpy, cools to form stars, 
especially in the innermost regions where the 
density is the highest. The radiation drag due 
to starlight acts on the cold gas, decreasing its 
angular momentum and causing an inflow into 
a reservoir around the central BH, to be sub- 
sequently accreted into it, increasing its mass 
and powering the nuclear activity. In turn, the 
feedbacks from SN explosions and from the 
active nucleus regulate the star formation rate 
and the gas inflow, and eventually unbind the 
residual gas, thus halting both the star forma- 
tion and the BH growth. Important parameters 
are the efficiency of SN energy transfer to the 



cold gas and the fraction of the QSO luminos- 
ity in winds. 

The model prescriptions are assumed to ap- 
ply to DM halos virializing at zvu ~ 1.5 and 

M vir > 2.5 x 10 n M G . These cuts are meant 
to crudely single out galactic halos associated 
with spheroidal galaxies. Disk (and irregular) 
galaxies are envisaged as associated primarily 
to halos virializing at z v ; r < 1 .5, some of which 

have incorporated most halos less massive than 
2.5 x 10 11 M virializing at earlier times, that 
may become the bulges of late type galaxies. 

The kinetic energy fed by supernovae is 
increasingly effective, with decreasing halo 
mass, in slowing down (and eventually halting) 
both the star formation and the gas accretion 
onto the central black hole. On the contrary, 
star formation and black hole growth proceed 
very effectively in the more massive halos, giv- 
ing rise to the bright SCUBA phase, until the 
energy injected by the active nucleus in the sur- 
rounding interstellar gas unbinds it, thus halt- 
ing both the star formation and the black hole 
growth (and establishing the observed relation- 
ship between black hole mass and stellar ve- 
locity dispersion or halo mass). Not only the 
black hole growth is faster in more massive ha- 
los, but also the feedback of the active nucleus 
on the interstellar medium is stronger, to the ef- 
fect of sweeping out such medium earlier, thus 
causing a shorter duration of the active star- 
formation phase. 

The basic yields of the model are the star- 
formation rate, tfr(t), as a function of the galac- 
tic age, f, (hence the evolution of the mass in 
stars, M sp h(f)X and the growth of the central 
BH mass, M Bll (t), for any given value of the 
halo mass, M v i r , and of the virialization red- 
shift, z v i r . These quantities are obtained solv- 
ing the system of differential equations given 
by Granato et al. (2004) and Cirasuolo et al. 

(2005) . 

5.2. Approximate analytic solutions 

A simplified set of equations, still keeping the 
basic physics, is presented by Shankar et al. 

(2006) , who derived an approximate analytic 
solution for the time evolution of the mass 
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Fig. 1. Comparison of the relationship between 
Mbh and Mh a i derived by Shankar et al. (2006, solid 
line; the thin horizontal bars show the uncertainties) 
with the prediction of the model by Granato et al. 
(2004; shaded area) and the observation-based esti- 
mates by Ferrarese (2002; dashed line) and Baes et 
al. (2003; dot-dashed line), referred to a virialization 
redshift z vil = 3. 

in stars, M stal , which, at the present time to, 
writes: 



M stal -(? ) = / su 



,/cosm M h (0) 

r 



(4) 



where / SU1V is the fraction of stars still sur- 
viving at the present time (/ sur v ~ 0.6 for a 
Salpeter IMF after about 10 Gyr from a burst), 
/cosm M h {0) is the initial gas mass (/ cosm ^0.17 
being the cosmic baryon to dark matter ratio 
and Mi, being the halo mass). Further, y = 
1 -R+a, R being the fraction of mass restituted 
by evolved stars (R » 0.3 for a Salpeter IMF), 
and 



Nsn fSN £sn 



(5) 



is effective efficiency for the removal of cold 
gas by the supernova feedback. Here /Vsn is 
the number of SNe per unit solar mass of con- 
densed stars, 6sn £sn the energy per SN used 
to remove the cold gas, and Eb the binding en- 
ergy of the gas within the DM halo, per unit 
gas mass. 

For z > 1 Shankar et al. (2006) obtain: 

a » 1.2(AW8xl0- 3 )(e SN /0.1) ■ 
• (£ SN /10 51 erg) [(l+z)^]- 1 • 



(M,/10 12 M Q ) 



-2/3 



(6) 



For large halo masses, where the stellar feed- 
back is less efficient (a ^ 1), the quantity 
1 - R + a is a slowly decreasing function of 
the halo mass, so that M stal - is approximately 
proportional to Mh- However, in this case, the 
fraction of gas turned into stars is controlled by 
the AGN feedback, which, as shown by the full 
treatment by Granato et al. (2004), for M h ^ 
3 x 10 11 M Q expels an approximately constant 
fraction of the initial gas, thus preserving the 
approximate proportionality between M star and 
Mh, in agreement with eq. Q. 

On the contrary, for M h <k 1O 12 M , y ^ 



a oc M 



-2/3 



and 



,/si 



,5/3 
l h ' 



(7) 



quantifying the effect of SN in decreasing the 
M star /M/, ratio with decreasing M/,. 

In the Granato et al. (2004) model, the ra- 
diation drag dissipates the angular momentum 
of the cool gas present in the central regions 
(Kawakatu & Umemura 2002). The gas then 
falls into a reservoir around the central BH at a 
rate 



M rt 



1.2 x 10"V(r)(i -O 



(8) 



where t is the effective optical depth of the cen- 
tral star forming regions. If most of the mass 
in the reservoir is ultimately accreted onto the 
central BH, and if other contributions to the 
BH mass can be neglected, we have: 



1.2xl0- 3 M star (l-O 



(9) 



Granato et al. (2004) assumed that the effective 
optical depth depends on the cold gas metallic- 



1 /3 

ity and mass t oc Z M„ as . The outcome of their 
numerical code yields, on average, Z oc M? 3 in 
the mass range 10 11 M e < M h < 3 X 10 13 M 
[cf. their Figs. (5) and (8)]. Since M gas ~ 
fcosmM h , one gets 



r oc M ] 



2/3 



(10) 



In the mass range Mf, ^ 10 11 M the optical 
depth is small (t <k 1); from eqs. and (II Ui 
we then obtain: 



M BH k Af star t oc M, 



7/3 



(ID 
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Thus this simple model predicts that the low 
mass slope of the Msa-Mf, relation is steeper 
than that of the M star - Mi, relation, because 
of the decrease of the optical depth with mass 

2/3 

t oc M h , entailing a lower capability of feed- 
ing the reservoir around the BH. 

For large halo masses r » 1 , and the model 
predicts Mbh k Mstar K Mi,. A roughly linear 
relationship between Mbh and M star (Mbh x 
M s ^ 2 6 ) has been observationally derived by 
Haring & Rix (2004). 

In Fig.[2we compare the Mbh-M/, relation 
derived by Shankar et al. (2006) with that of 
Ferrarese (2002), who first investigated this is- 
sue from an observational point of view. She 
derived a power-law relationship between the 
bulge velocity dispersion and the circular ve- 
locity, v c , for a sample of spiral and elliptical 
galaxies spanning the range 100 5= v c ^ 300 
km s , and combined it with an approximated 
relationship between v c and the virial veloc- 
ity, v v i r , based on the numerical simulations by 
Bullock et al. (2001) and with the M/,-v vir rela- 
tionship given by the ACDM model of the lat- 
ter authors for a virialization redshift z v i r ~ 0, 
to obtain a Mi-cr relation. Coupling it with one 
version of the observed BH mass vs. stellar ve- 
locity dispersion relationship (Mbh k cr 4 58 ) 
she obtained M B h oc M" x , with a = 1.65-1.82. 
Baes et al. (2003) with the same method, but 
assuming Mbh 00 tr 4 02 and with new veloc- 
ity dispersion measurements of spiral galax- 
ies with extended rotation curves, yielding a 
slightly different v c -cr relation, found Mbh k 
Ml- 27 . 

As shown by Fig. 6a of Shankar et al. 
(2006), their relationships lie substantially be- 
low ours. 1 It should be noted, however, that 
the M/,-v v i r relation depends on the virializa- 
tion redshift. For z v ; r ^ 3, its coefficient would 
be a factor of ^ 4.25 lower than that used by 
Ferrarese (2002) and Baes et al. (2003) and the 
coefficients of the Mbh 00 M/, relations would 
be larger by a factor of 5.6 in the case of 



1 Both Ferrarese (2002) and Baes et al. (2003) 
adopted V m cc [n M (A vlr /200)]°' with a = 1/3, fol- 
lowing Bullock et al. (2001), while the correct value 
is a = 1/2. This correction has, however, a minor 
effect. 



eq. (6) of Ferrarese(2002) or of ^ 4 in the 
case of Baes et al. (2003), bringing them much 
closer to the results by Shankar et al. (2006; see 
Fig.Q, and to the predictions of the Granato et 
al. (2004) model. 

5.3. The evolution of nuclear activity 

According to the Granato et al. (2004) model 
the SMBH growth parallels the star forma- 
tion in spheroidal galaxies. However, the BH 
growth is much faster than the increase of the 
mass in stars. In fact, the timescale for the 
former is roughly the e-folding timescale for 
Eddington-limited accretion, f e f = fg e/( 1 — e) 
(where is the Eddington time and e is the 
BH mass to energy conversion efficiency). For 
e = 0.1, f e f - 5 x 10 7 yr. The star forma- 
tion timescale, which is the shorter between the 
cooling and the free-fall times, increases from 
the denser inner regions outwards (where most 
of the mass is); its effective values are typi- 
cally of a few to several hundred million years. 
This means that the MBH/M stal - ratio increases 
rapidly until the AGN feedback sweeps out the 
residual gas. 

The exponential growth of the active nu- 
cleus implies that its bolometric luminosity is 
very low except in the final several e-folding 
times. Also the nucleus is surrounded by a 
dusty environment with a large column den- 
sity, even if we ignore the contribution of the 
obscuring torus. For example, if we adopt for 
the gas distribution a King profile with core ra- 
dius ~ 0.3-1 kpc, we obtain, for galaxies with 
Mi, i 2 1O 12 M , column densities to the nu- 
cleus of N H ~ a few x 1 23 to a few x 1 24 cm" 2 
(Granato et al. 2006). Correspondingly, we ex- 
pect large optical extinctions Ay ~ 200(/V///3x 
10 23 cm" 2 ). 

The nuclear activity then shows up first in 
hard X-rays. As noted in §|4] the X-ray visibil- 
ity times are in the range 0.1-0.3 Gyr, corre- 
sponding to the last 2-6 e-folding times of the 
BH growth. In this period, the M B H/M star ra- 
tio can be <k 1 . The AGN feedback will then 
progressively sweep out the residual gas, de- 
creasing the extinction until the AGN enters its 
optical visibility phase, whose duration is ex- 
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Fig. 2. Contribution of SMBH associated to spheroidal galaxies at z — 2 (solid line) to the local SMBH 
mass function (dashed line; the shaded area shows the uncertainty) estimated by Shankar et al. (2004), 
assuming a M BH /M stilr ratio 3 times higher than in the local universe and adopting a M/L ratio appropriate 
for a Salpeter (left-hand panel) IMF or an IMF flattening below 1 M Q (right-hand panel), with a dispersion 
of0.3dex. 



pected to be a substantial fraction of the last 
e-folding time. 

These expectations are fully borne out by 
the study of 20 SCUBA galaxies brighter than 
4 mJy at 850 /mi, lying in the region of the 
2Ms CDF-N Chandra observations (Alexander 
et al. 2003, 2005a,b). As shown by Granato 
et al. (2006), the data imply accretion rates 
and column densities fully compatible with 
the predictions of the Granato et al. (2004) 
model. Borys et al. (2005) estimate, for these 
objects, Mbm/M^m ratios one or two orders 
of magnitude smaller than local galaxies of 
comparable stellar mass, suggesting a subse- 
quent growth lasting 2 to 5 e-folding times, at 
the Eddington limit. The lack of objects with 
M BH /M star > 0.1(M BH /M ster )i ocal may suggest 

that in the last ~ 2-3 e-folding times the AGN 
feedback begins quenching the star formation 
rate so that the sub-mm emission decreases and 
galaxies slip below the 850 //m flux limit. 

A different history of the coevolution of 
SMBHs and host galaxies is apparently im- 
plied by the study by Peng et al. (2006), who 
estimated the Mgu/M^ ratios for a sample of 
optical quasars at z ~ 2 and found them to 
be 3-6 times larger than in the local universe. 
This conclusion however relies on a number 
of very uncertain assumptions, as discussed by 
the authors themselves. A further constraint 
come from the SMBH local mass function (see 



§[4}. An estimate of the galaxy stellar mass 
function at z - 2 has been obtained by Drory 
et al. (2005). Others can be obtained from the 
galaxy luminosity functions at a similar red- 
shift derived by various groups (Saracco et al. 
2006; Dahlen et al. 2005; Fontana et al. 2004; 
Caputi et al. 2005). As shown by Fig. |2j a 
MBn/M stm - ratio 3 times higher than in the lo- 
cal universe would imply that, at large SMBH 
masses, there is no room for growth at z £ 2. 

The constraint is attenuated if the IMF flattens 
to a slope of -0.4 below 1 M e (see Romano et 
al. 2002 and Granato et al. 2004), decreasing 
the M/L ratio in the ZT-band to ^ 0.3 (this ratio 
is =! 0.5 for a Salpeter IMF). Nevertheless, even 
in this case a rather large contribution to the lo- 
cal SMBH mass function is already in place at 
z - 2, and this is very difficult to reconcile with 
the optical and X-ray data on the evolution of 
the AGN luminosity function. 

6. Conclusions 

Huge amounts of data have been accumulating 
in recent years, providing rather detailed in- 
formation on the mass-dependent evolution of 
galaxies and of active nuclei. Rather unexpect- 
edly, it turned out that intense star formation 
shifts from massive galaxies at high redshifts 
to lower mass galaxies at later time, a pat- 
tern referred to as "downsizing" (Cowie et al. 



De Zotti et al.: Co-evolution of QSOs and spheroidal galaxies 



289 



1996). Granato et al. (2001, 2004) showed that 
this evolutionary behaviour, which apparently 
contradicts expectations from the hierarchical 
structure formation paradigm, can be recon- 
ciled with it if the main driver of evolution is 
not the merging sequence of dark matter ha- 
los, but are baryon processes, involving strong 
feedback effects from star formation (and su- 
pernova explosions) and from nuclear activ- 
ity (Anti-hierarchical Baryon Collapse - ABC 
- scenario). 

An implication of this scenario is that the 
luminosity function of massive galaxies under- 
went essentially passive evolution since z — 
1.5, little affected even by the "dry mergers" 
advocated by several authors (e.g. Bell et al. 
2005; Faber et al. 2006; van Dokkum 2005), 
but strongly constrained by the DEEP2 Galaxy 
Redshift Survey (Bundy et al. 2006). The ABC 
scenario also explains a broad variety of ob- 
servational data, as shown by Granato et al. 

(2004) , Silva et al. (2005), Cirasuolo et al. 

(2005) , and Shankar et al. (2006). We have re- 
called here the close match of the virial veloc- 
ity function with the velocity dispersion func- 
tion of galaxies, and the M B h-o~ relation, both 
determined from data on local galaxies, but 
whose physical explanations trace their origin 
to substantial redshifts. 

In this framework, also the AGN evolu- 
tion is tightly connected with star formation, 
rather than with mergers (although star forma- 
tion may be triggered by merging events), and 
the BH growth is dominated by radiative ac- 
cretion, rather than, e.g., by coalescence of BH 
associated to merging sub-units. This implies 
that during the intense star-formation phases, 
characterizing sub-mm bright galaxies, the BH 
masses are still below, even by large factors, 
the present day values, consistent with the find- 
ings by Alexander et al. (2003, 2005a,b) and 
Borys et al. (2005). Also, the detectable QSO 
activity is delayed compared to the onset of 
vigorous star formation, since: a) the BH mass 
and its bolometric luminosity grow exponen- 
tially on a timescale f e f ^ 5 x 10 7 yr and are 
therefore small except in the final few tenths of 
Gyr before reaching the final mass; b) the nu- 
clear activity is highly obscured even in hard 



X-rays before the AGN feedback begins to 
push away the residual dusty medium. 
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